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AN ANALYSIS OF LAMINAR FREE - CON VE CTION FLOW AND HEAT TRANSFER 
ABOUT A FLAT PLATE PARALLEL TO THE DIRECTION 
OF THE GENERATING BODY FORCE 1 

By Simon Ostbach 


SUMMARY 

The free-convection flow and heat transfer (generated by a 
body force) about aflat plate parallel to the direction of the body 
force are formally analyzed and the type of flow is found to be 
dependent on the Orashof number alone. For large Orashof 
numbers (which are of interest in aeronautics), the flow is of 
the boundary-layer type and the problem is reduced in a formal 
manner, which is analogous to Prandtl’s forced-flow boundary- 
layer theory, to the simultaneous solution of two ordinary 
differential equations subject to the proper boundary conditions. 

Velocity and temperature distributions for Prandtl numbers 
of 0.01, 0.72, 0.7S8, 1, 2, 10, 100, and 1000 are computed, and 
it is shown that velocities and Nusselt numbers of the order of 
magnitude of those encountered in forced-convection flows may 
be obtained in free-convection flows. The theoretical and 
experimental velocity and temperature distributions are in good 
agreement. 

A flow and a heat-transfer parameter, from which the impor- 
tant physical quantities such as shear stress and heat-transfer 
rate can be computed, are derived as functions of Prandtl num- 
ber alone. Comparison of theoretically computed values of the 
heat-transfer parameter with values obtained from an approxi- 
mate calculation and experiments yielded good agreement over a 
large range of Prandtl number. Agreement between the 
theoretical values and those obtained from a frequently used 
semiempirical heat-transfer law was good only in restricted 
Prandtl number ranges ( depending on an arbitrary constant) . 

INTRODUCTION 

Two important types of fluid flow problems involving heat 
transfer are those of forced and those of free convection. 
By forced-convection flow is meant flows maintained mechan- 
ically as, for example, by a pressure drop or ah agitator. 
Free-convection flow, on the other hand, results from the 
action of body forces on the fluid, that is, forces which are 
proportional to the mass or the density of the fluid. The 
flow is generally produced in the following manner: Consider, 
for example, a fixed object (such as a plate) in a quiescent 
fluid subject to a body force. When the plate is at the same 
temperature as the surrounding fluid, the body forces acting 
on the fluid are in equilibrium with the hydrostatic pressure 
and no flow ensues in the steady state. If a temperature 
gradient normal to the body force is imposed by heating 
(or cooling) the plate, there will exist a defect (or excess) of 


body force because of the decreased (or increased) density, 
with the fluid closer to the plate having the greater defect 
(or excess) than that away from the plate. This unbalance 
of the forces causes the fluid to be accelerated with the 
particles nearer the plate moving more rapidly than those 
farther from the plate. Free-convection flow has usually 
been considered to be generated in a gravitational field where 
the previously mentioned defect or excess of body force was 
the Archemedian (buoyancy) foice. However, since cen- 
trifugal forces are also proportional to the fluid density, free- 
convection flows can also be set up by the action of such 
forces. (See ref. 1 for a more explicit discussion of the 
development of. free-convection- flows by centrifugal forces.) 

Free-convection flows produced by centrifugal forces are 
now of practical importance in aeronautics because many 
aircraft propulsion systems contain components (such as gas 
turbines and helicopter ram jets) which rotate at high 
speeds and in which heat is being transferred. The method 
of free-convection cooling of gas-turbine rotor blades where 
the centrifugal forces create a free-convection flow of the 
coolant in the blade passages is an example of a practical 
application of the free-convection phenomenon in aero- 
nautics. Also, free-convection flow due to centrifugal force 
is superimposed on the flow through helicopter ram jets and 
on the flow of cooling air in hollow rotor blades of air-cooled 
turbines and, under proper conditions, can appreciably 
influence the resultant flow and heat transfer. 

• Ah a simplification of the many free-convection problems 
which are now of some consequence in aeronautics, consid- 
eration is here given to the special case of free-convection 
flow about a flat plate parallel to the direction of the gener- 
ating body force. The experimental and theoretical consid- 
erations of Schmidt and Beckmann (ref. 2) concerning the 
free-convection flow of air subject to the gravitational force 
about a vertical flat plate constitute the most complete 
treatment of this subject up to the present time. Eckert 
(ref. 1) as well as others has further verified and extended 
the experimental results of Schmidt and Beckmann, and 
Schuh (ref. 3) has extended the numerical calculations by 
computing the velocity and temperature distributions for 
several Prandtl numbers different from that for air. How- 
ever, all the theoretical work in these references is based on 
the incompressible equations in which the density (or tem- 
perature) variation is introduced in the buoyancy term alone. 


i Supersedes NAOA TN 2635, "An Analysis of T-nmlnar Free- Convection Flow and Heat Transfer about a Flat Plate Parallel to the Direction of the Generating Body Force*' by Simon 
Oatracb, 1052. 
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Various terms are omitted from the equations at the start 
on the basis of either intuitive arguments or no arguments 
at all. Although a theoretical development made in such a 
manner led to good final results, the significance of all the 
important factors associated with the free-convection flow 
phenomenon is not obtained from such an analysis. 

The problem of free-convection flow as produced by a body 
force about a flat plate in the direction of the body force was 
studied at the NACA Lewis laboratory during 1951 and is 
treated in a formal and more general manner herein. The 
method used is somewhat similar to that used in reference 4 
wherein consideration was given to the free-convection flow 
at hig h Grashof numbers in a horizontal cylinder which had a 
variable surface-temperature distribution. The application 
of this method to the present problem leads to a development 
which is analogous to Prandtl’s treatment of high Reynolds 
number forced-convection flows. Although the final equa- 
tions obtained by this method are the same as those of 
Schmidt and Beckmann, this more general approach not only 
clearly demonstrates the significance of all the important 
parameters and assumptions and hence leads to a better 
understanding of this type of flow but also indicates the 
quantitative limitations of the theory. In addition, the 
numerical solutions of references 2 and 3 are herein extended 
to cover a more complete range of parameters. The new 
calculations yield information on the free-convection flow for 
Prandtl numbers corresponding to those of liquid metals, 
gases, liquids, and very viscous fluids. 

ANALYSIS 

STATEMENT OF PROBLEM AND BASIC EQUATIONS 

The steady-state equations expressing the conservation of 
mass, momentum, and energy for a compressible, viscous, 
and heat-conducting fluid subject to a body force together 
with an equation of state govern the flow and associated 
temperature distribution aboyt the plate. These equations 
in Cartesian tensor notation are (see ref. 5), respectively, 
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(A complete list of the symbols used herein is given in 
appendix A.) Por the two-dimensional case, equations (1) 
to (4) represent a system of five equations in the five depend- 
ent variables Hi, U 2> p, P, and T. Por later use, equation (4) 
can be written 


where K. and 0 are the coefficients of isothermal compressi- 
bility and volumetric expansion, respectively (see ref. 6). In 
addition to a general state equation, such as is given in 
equations (4) or (4a), it will be convenient at timeB in the 
discussion to refer to some specific state equation. To this 
end, the'equation of state for an ideal gas 

P—pUT (4b) 

will be used. 

Particular consideration is here given to the two-dimen- 
sional free-convection flow about a semi-infinite vortical flat 
plate. The A^axis of the coordinate system is taken along 
the plate and the Aj-axis, normal to it. No distinction is 
made as to the specific type of body force acting, for example, 
gravitational or centrifugal, but the force is assumed to be 
acting in the vertical direction only (that is, parallel to the 
plate). Centrifugal and Coriolis forces which are connected 
with flows on curved paths and with rotating systems 
generally vary with position and velocity. However, in 
order not to make the analysis unduly complicated, the body 
force is taken to be constant. 

In order to define the problem clearly, a choice must still 
be made of the position of the origin of the coordinate sys- 
tem. Before making a definite decision on this point, note 
that for constant plate temperatures there are four permu- 
tations of the body-force direction (either upward or down- 
ward) and the plate thermal condition (either heated or 
cooled) which will lead to free-convection flows. Once the 
position of the edge of the plate, which is also to be the 
origin of the coordinate system, is decided, there are two 
combinations of the body-force direction and plate thermal 
condition that will yield flows which proceed away from the 
edge. It is this type of flow that is amenable to the type of 
analysis to be made here. This point will be more fully disi 
cussed subsequently. If the edge of the plate (recall that a 
semi-infinite plate has but one edge) is taken at the bottom of 
the plate (that is, the plate extends to + <» in the Aj-direc- 
tion), the two combinations loading to flows in the proper 
direction (upward in this case) are, respectively, the body 
force acting downward with a heated plate and the body 
force acting upward with a cooled plate. The equations 
developed for one of the cases reduce directly to those for the 
other. The remaining two permutations, namely, the body 
force acting downward with a cooled plate and the body force 
acting upward with a heated plate, would yield flows which 
proceed downward or toward the edge of the plate if this 
edge were taken at the bottom of the plate. This type of 
flow would violate a physical condition of the problem which 
states that the flow starts at the plate edge. The latter 
combinations hence will not be considered further. 

Because the two acceptable configurations can be reduced 
essentially to one, for the development to be given here, the 
origin of the coordinate system will be taken at the bottom 
of a heated plate, with the body force acting downward. 
The assumption is now made that the viscosity and thormal- 
conductivity coefficients are functions of the temperature' 
alone and obey the following laws: 


dp=p(A dP — $ dT) 


(4a) 
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variation in tlie fluid, tlie free-convection velocity should 
then depend directly on e, and the variations in pressure and 
density (from the static, e=0, case) due < to the temperature 
differences should also depend on 6. Thus 


The choice of the body-force direction together with equa- 
tions (5) alters equations (2) and (3) so that they become 
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Note that the only nonzero component of the body force is 
the Xi-component. 

BOUNDARY CONDITIONS 

The boundary conditions associated with the given prob- 
lem are that: 

(a) The fluid must adhere to the plate (the no-slip condi- 
tion of viscous flows) and the plate must be a streamline, 
or mathematically, 


C7i(-^"i>0)= C7j(Xi,0)=0 (8) 


(b) The temperature of the fluid at the plate must he equal 
to the plate temperature, that is, 


T{X h 0)=T, (9) 


(c) The velocity TJ\ at large distances from the plate must 
be undisturbed, or 


Ux(X u »)=0 


( 10 ) 


(d) The temperature at large distances from the plate 
must be equal to the undisturbed fluid temperature, or 


TO,»)=jp. 

SIMPLIFICATION OF EQUATIONS 

Lot a small quantity c now be defined as 


( 11 ) 
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where —f x denotes the i?) -component of the body force per 
unit mass, u {> a, <p, and 6 denote dimensionless functions 
(which, in general, can be functions of e), l is some character- 
istic length (for example, the distance from the edge of the 
plate to the point of interest), P, and p, are the pressure and 
the density, respectively, for the static case (77^0 or 
c — 0), and P a and p„ denote constant values of the pressure 
and the density (that is, the values if no force field were 
present) defined by the state equation (in the case of a gas, 
in particular, P a ,=p a ,RTco)‘. Because there is no character- 
istic velocity associated with the type of flow under con- 
sideration, the velocity is dimensionahzed by the factor 
given in parentheses on the right side of equation (13). 

In order to determine the static quantities, it will at first 
be convenient to consider the particular case of a gas. The 
problem is then considered with the temperature uniform 
throughout the flow field at the value (therefore there, 
will be no flow and 77<=0). For this situation, equations 
(4b) and (6) become 
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(It should be noted that eq. (18) expresses the physical fact 
previously stated that the body force and hydrostatic 
pressure are in equilibrium for the static case.) Substitution 
of equation (17) into equation (18) leads to 

P.=P„ex p(-^^i) (19) 

and equation (19) together with equation (17) and the 
equation defining P m and p„ yields 


which is a measure of the magnitude of temperature variation 
in the flow field. The coefficient of volumetric expansion j3 
is generally of the order of magnitude between 10 -2 and 10 -4 
(see table 15 of ref. 7, for example) and for gases, p=l/T. 
(Thus, for gases, if (3 is taken to he constant, e= (T Q — T„)jT a ; 
that is, e is the relative temperature difference.) The 
coefficient /3 will be assumed constant. Because in the 
Bteady state flow ensues only when there is a temperature 


p,=p„exp(— ^-Xi^=p„exp^— Xj) (20) 

If the exponential in equation (20) is expressed in terms of 
its series expansion, that equation becomes 

p,= p„( 1-^^+...) (21) 
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A computation of the second term in the parentheses of 
equation (21) for the case of air under normal conditions 
•with fx—g and the fact that X x is of the order of magnitude 
l show that p w gllP a ~ 10~ s Z/foot. For the type of problem 
under consideration, l ■will always be of unit order of magni- 
tude so that even if the body force f x represents a cen- 
trifugal force many times that of gravity, the inequality 
p»/a-Xi/P„«: 1 may stillbe satisfied. Thus, in the subsequent 
development it will be assumed that This assump- 

tion, which was justified by the computation for the case of 
a gas, is expected to be reasonable for other fluids as well. 
The physical interpretation of this assumption is that under 
static conditions (e=0), the density (or pressure) is not 
affected by the force field. 

In order that all quantities in the following equations be 
dimensionless, it is further necessary to define x { =X t fl, 
where the x t are now dimensionless space coordinates. Sub- 
stituting these new coordinates along with equations (13) to 
■ (16) into equations (1), (6), (7), and (4a) and noting equations 
(18) and that yield, on neglection of terms of higher 

order in e compared with those of order e, 
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Thus, to a first approximation, equations (22) to (26) to- 
gether with the boundary conditions replace the original 
equations and boundary conditions. (Note that for gases, 
j8T„=l.) 

The prime assumption made in this analysis is that the 
higher-order terms in e are negligible, which implies that e 
is small, and consequently, that the temperature difference 
or 0 is moderately small. It is a consequence of this assump- 
tion alone that the basic equations were simplified to equa- 
tions (22) to (26)- wherein tdie viscosity term in the energy 
equation is neglected and the only coupling of the momentum 
and energy equations occurs by means of the body-force 
term in equation (23). As a result of this assumption, the 
variations of the viscosity and heat-conductivity coefficients 
with temperature are also negligible. Without any discus- 
sion, the authors of reference 2 start directly from simplified 
equations of the same form wherein the pressure terms in the 
energy equation were also neglected. In reference 3 some 
intuitive arguments are given to justify the simplified 
equations. 

It is now convenient to revert to the more familiar notation 
where x=x u y=x 2 , u—u u and ti=u%. Equation (22) implies 
the existence of a stream function ^ such that 



dip—KP^da — pT^dO ( 26 ^ 

where P»/p»/x^=W 6V and the Grashof number Gr and the 
Prandtl number Pr are defined as 


and 


Gr= 





Physically speaking, the Grashof number represents the 
ratio of the body forces to the viscous forces. The free- 
convection flows of interest here are those associated with 
large Grashof numbers. The factors K and ($ in equation 
(26) may well be taken to be constants (see ref. 6). 

The boundary conditions (eqs. (8) to (11)) in terms of the 
new dimensionless variables are 
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where subscripts denote differentiation. Applying equation 


(31) to equations (23) to (25) yields, respectively, 

~ <p)=4> v \p xr —ip I \l/ n ,+Ncr x (32) 

— ^( A ^*) = ~ Mrc+'f'x'l'xt+Nffp (33) 

^p^(A (34) 

The boundary conditions (eqs. (27) to (30)) become 

'f'v (*»0)= ife(*,0)= 0 (36) 

e(x ’ 0)= w: (36) 

t v (x,co)=e{x,<*)=0 (37) 


Equations (32) to (34) and equation (26) form the system 
of equations for the four unknown functions 8, <p, and a 
of the problem. The system is nonlinear, and therefore 
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further simplification of the equations would he desirable. 
Just as in the case of forced-convection flows where the 
Reynolds number determines the type of flow or, in mathe- 
matical terms, the type of solution, the Grashof number is 
the prime factor for free-convection flows. For the case of 
small Grashof number, it can be seen from equations (32) 
to (34) that a perturbation in the small parameter Or will 
yield a system of linear equations. For Grashof numbers 
of unit order of magnitude, no further important simplifica- 
tion can be made and the solutions would have to be obtained 
numerically. For the other limiting case, that of large 
Grashof numbers (which is- the case under consideration 
herein), it would, at first thought, appear that some simpli- 
fication could be obtained by performing a perturbation in 
the small parameter 1/Gr. However, this would then 
imply that the term containing the highest-order derivatives 
(the left term in equations (32) to (34)) could, among others, 
be neglected. (This argument would also imply that the 
body-force term <p in equation (32), which is essentially, 
causing the flow, could also be neglected.) The omission 
of the highest-order derivatives from consideration, however, 
would lead to solutions which would not satisfy all the 
boundary conditions. Problems of this type are referred 
to as singular perturbation problems. For further dis- 
cussions of singular perturbation problems, see references 
8 and 9. 

Equations in which a small parameter multiplies the 
highest-order terms are said to be of the boundary-layer 
type, because in order for solutions which satisfy all the 
boundary conditions to be obtained, the highest-order terms 
must be considered near the boundary. This fact implies 
the existence of a thin region, called the boundary layer, 
wherein the functions vary rapidly from the value at the 
boundary to that in the flow outside this layer. The con- 
clusion to be drawn from the preceding discussion is that 
for large Grashof numbers the flow is of the boundary- 
layer type. Schmidt and Beckmann (ref. 2) also made 
the boundary-layer assumptions in their theoretical develop- 
ment, and these assumptions were justified on the basis of 
their experimental observations. The Grashof numbers for 
their experiments were of the order of 8X10 8 . 

In view of the fact, previously discussed, that highest- 
order derivatives of each dependent variable as well as of 
those terms of physical importance (as, for example, the 
body-force term) must be retained in the boundary layer, it 
is convenient to make both sides of each of the equations 
of the same order in Gr. In this way, as will be shown, the 
equations will be further simplified. It is thus convenient 
to make the following transformations in the system of 
equations (32) to (34) and (26) and then to retain only the 
dominant parts (that is, those multiplied by Gr to the 
highest power) of each individual term. 

Let y—Gr T y, a=Gr'a, ¥>=$5, and 8=6. Then 

equations (32) to (34) and (26) become 
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It now can be seen that by proper choice of r, s, and t a 
transformation of the type given provides a means for 
making the important terms in the differential equations of 
the same order in Gr. Thus if r=%, 8=—%, and t=— 1, 
equations (38) to (41) become 

fe— <p— fete+NoF, 

(42) 

Nr;= 0 ( Gr 5 )e0 

(43) 


(44) 

d<o-\- pT (O cZ0=O 
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More generally, if N is very much different from unit order 
of magnitude, a value of t can always be chosen (depending 
on N) such that equations (42) to (45) are obtained. (For 
any negative t less than —1, the last term of eq. (42) 
will also disappear.) 

There are now several important points to be discussed 
concerning the transformation just made and the resulting 
simplified equations. First, it should be noted that the 
transformation is merely a formal expression of the bound- 
ary-layer assumptions first made by Prandtl and hence the 
solutions will be asymptotic for largo Gr. Second, the 
second equation of motion here also reduces to state that 
the pressure across the boundary layer is constant. Third, 
the pressure terms in the energy and state equations are 
here found to be negligible. This fact verifies a priori 
assumptions made by others from the physics of the problem. 
Finally, note that integration of the general state equation 
(independent of pressure) as now given by equation (45) 
leads to 

V+PTJ=0 (46) 

where the constant of integration has been taken as zero 
without any loss of generality. For the particular case of 
a gas, fi=l /T a so that equation (46) becomes 

<p- (- 0=0 

The boundary conditions (eqs. (35) to (37)) now can be 
written 


i^(r,0)=^(x,0)=0 

(47) 

0(z,°)=^ 

(48) 

^(x,^)=e{x,co)=o 

(49) 


If now it is assumed that <F*=0 in equation (42) since con- 
sideration is here being given to a flat plate, and if Ip is 
eliminated from equation (42) by use of equation (46), 
there results the system of equations 
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'Pzz+ffT„6=ifolr x ;—Tp I ifa (50) 

/ 
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Thus the problem has been reduced to the solution of the 
two simultaneous partial differential equations (eqs. (50) and 
(51)) subject to the boundary conditions (eqs. (47) to (49)). 

Final simplification of the equations is made by applica- 
tion of the so-called similarity transformation of boundary- 
layer theory. Thus, let 


and 



^=(4x)*i ? (ij) 


(52) 

(53) 



(54) 


Then equations (50) and (51) are reduced to the following 
ordinary differential equations : 

F'"+3FF"-2F' 2 +H=0 (55) 


Schmidt-Beckmann calculations wherein the value of Prandtl 
number of 0.733 was used. The particular values of the 
Prandtl numbers given were chosen to correspond to those 
for liquid metals, gases, liquids (such as water and oil), and 
very viscous liquids (such as glycerin or oils at very low 
temperatures). 

RESULTS 

VELOCITY AND TEMPERATURE DISTRIBUTIONS 


By means of the various transformations made in the 
analysis it can easily be verified that 
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where 
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(62) 


H"+2PrFH'=0 (56) 

where the primes denote differentiations with respect to 17 . 
The reciprocal one-fourth power similarity as given in equa- 
tion (52) is characteristic of free-convection flows just as the 
reciprocal square-root type is characteristic of the forced- 
convection flows. The boundary conditions become 

F'(0)=F(0)=0 (57) 

17(0) = 1 (58) 

F'(oo)=ff («)= 0 (59) 

The use of a transformation like equation (52) essentially 
specifies an additional boundary condition, namely, that the 
conditions to be satisfied at y= °o (or 17 = <») should also be 
satisfied at x= 0 . It is for this reason that the flows pre- 
viously discussed which would flow toward the edge (down- 
ward) are not amenable to this type of analysis, for such flows 
would violate this additional condition, which essentially 
states that the boundary-layer development starts at the 
edge of the plate. 

SOLUTION OF THE BOUNDARY- VALUE PROBLEM 

The solutions of the simplified equations (55) and (56), ’ 
satisfying the boundary conditions as given by equations 
(57) to (59), were obtained by use of an IBM Card-Pro- 
grammed Electronic Calculator. A detailed account of the 
procedure followed in the determination of the unknown 
functions is presented in appendix B by Dr. Lynn U. Albers. 
The functions F and H together with their derivatives are 
given in table I for Prandtl numbers of 0 . 01 , 0.72, 0.733, 1 , 2 , 
10, 100, and 1000. Even though the Prandtl number for air 
is taken as 0.72 in this report, the solutions for Pr=0.733 
were also computed and are presented as a check with the 


Equations (60) to (62) relate the physical quantities to the 
dimensionless functions F and H which are now known. 
The dimensionless velocity and temperature distributions as 
given by equations (60) and (61) are presented in figures 1 
and 2 , respectively, as functions of 17 for the various values 
of Prandtl number. The computations made here agree 
with those for Pr~ 0.733 as given in reference 2 up to the 
third significant figure. For Pr— 10 , 100 , and 1000 , the 
present results agree in general with those of reference 3. 
Since only curves are presented in reference 3, the precision, 
of the agreement cannot be stated. 

The maximum values of the dimonsionless velocity dis- 
tributions occur at larger values of the argument 17 as the 
Prandtl number decreases and the velocities decrease with 
increasing Pr. It should also be noted that tho dynamic 
and thermal boundary-layer thicknesses can he estimated 
from the abscissas of figures 1 and 2 , respectively, and that 
for Pr> 1 the velocity boundary layer is much thicker than 
the thermal boundary layer. 

The occurrence of /x (or Qr x as given by eq. (60)), which 
may be very large for flows generated by cen trifugal forces, 
in the denominator of the ordinate implies that velocities of 
appreciable magnitude can be associated with suoh free- 
convection flows. In particular, if J x = 10 8 feet per second 
squared, which is a reasonable conservative figure for pres- 
ent-day rotating systems, e= 0.2 (which is within the limits 
cf the theory presented herein), and arbitrarily X=0.25 
foot, then the maximum velocity attained at a Prandtl 
number of 0.72 is approximately 125 feet per second. This 
value of the maximum velocity could, of course, be doubled 
or even tripled under the proper condi tions. One limitation 
to a calculation of this sort, as can be seen by comparison of 
the denominators of the left and middle terms of equation 
(60), should be kept in mind; namely, the limiting Grashof 
number for laminar flows. In lieu of a complete stability 
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analysis ort tliis type of flow, this limiting value is taken to be 
10 9 , as indicated in reference 10. Consideration of this lim- 
itation then .implies (see eq. (60)) that for large laminar 
velocities either v a must be large or X must be small. 

COMPARISON WITH EXPERIMENTS 

Careful experiments of free-convection flows (as generated 
by gravitational forces) about vertical flat plates were made 
by Schmidt and Beckmann (ref. 2) in which velocity measure- 
ments at various points along the plate were made by means 
of a quartz-filament anemometer and the temperature meas- 
urements were obtained by means of manganese-constantan 
thermocouples. Eckert (ref. 1) performed similar experi- 
ments in which the measurements were made by means of a 



Zehnder-Mach interferometer. The results of both sets of 
experiments are in good agreement, but since the data pre- 
sented in reference 2 by Schmidt and Beckmann appear in 
more detail, these data will be used for comparison with the 
theory. • 

The experiments of reference 2 were performed on two 
different (in that the edges were smoothed either symmetri- 
cally or not) 12- by 25-centimeter plates and on one 50- by 
50-centimeter plate. It should here be pointed out that the 
results for the two smaller plates were almost identical and 
that the flow was entirely laminar except near the outer 
edge of the boundary layer where the slight turbulence of 
the room air disturbed the measurements somewhat. (This 
effect was also observed by Eckert.) Large periodic oscilla- 
tions of the flow near the downstream edge of the larger plate 
were observed in addition to the slight turbulence near the 



Figure 1 .—Dimensionless velocity distributions for various Prandtl numbers. 


Figure 2. — Dimensionless temperature distributions tor various Prandtl numbers. 
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outer edge of the boundary layer. Hence the data from the 
larger plate should not he expected to yield completely satis- 
factory agreement with the laminar theory as presented 
here. 

Since the physical quantities can be expressed in terms of 
a single variable as in equations (60) and (61), it is to be ex- 
pected that the data taken at the various points along the 
plates should all lie on a single line if the data are correlated 
according to equations (60) and (61). Thus for the smaller 
plates where (T 0 -T a> )=95.22° R and T a =518.68° R, equa- 
tions (60) to (62) become 



(63) 

T- .518.68 , 

95.22 ~ HM 

(64) 

nn Y 


77 = 88.26-=^ 

X i 

(65) 


The velocity and temperature distributions are so plotted in 
figures 3 and. 4, respectively, as are the curves computed 
theoretically for Pr= 0.72. It can be seen that the agreement 
is in general very good for small values of 17 and somewhat 
less satisfactory though still rather good for the larger values 
of 17 . The scatter in the range of the larger values of 1 ? is 
believed to be caused by the previously discussed room 



Ftqoee 3.— Comparison of small plate experimental and theoretical velocity distributions 
for Prandtl ntnnber of 0.72. 



Figure 4. — Comparison of small plate experimental and theoretical tomporaturo distribu- 
tions for FrandU number of 0.72. 


turbulence. It should also be noted that the points farthest 
away from the theoretical are those measured near the 
leading edge. These points should not, of course, be expected 
to agree- too well with the theory since the boundary-layer 
assumptions made in the theoretical development imply that 
the distance along the plate is large as compared with the 
boundary-layer thickness. Hence, this assumption is invalid 
near the leading edge. Schmidt and Beckmann obtained 
closer agreement between the theory and the experiments for 
the temperature data and poorer agreement for the velocity 
data by basing the kinematic viscosity coefficient in equation 
(62) on the plate temperature rather than on' the undisturbed 
stream temperature as was done here. 

Tor the larger plate, (T 0 — T») = 83.7° R and T a =627.14° 
R, so that equations (60) to (62) become 


4.522-v^ F ^ 

( 66 ) 

T~ 527.14 
83.7 

(67) 

77=83.93 

( 68 ) 



The velocity and temperature distributions for this experi- 
ment are plotted in figures 5 and 6 , respectively, and again 
the theoretical curves for Pr=0.72 are included. In figure 5 it 
can be seen that for large 17 the agreement is rather poor, 
particularly for the data for both small and large values of X. 
The poor agreement for small values of X is again due to the 
theory limitation near the edge of the plate and for large 
values of X, to the fact that the flow was becoming turbulent 
'there. 





FREE-CONVECTION FLOW AND HEAT TRANSFER ABOUT A FLAT PLATE 


71 


FLOW AND HEAT-TRANSFEE PARAMETERS 

In addition to the velocity and temperature distributions, 
it is often desirable to compute other physically important 
quantities (such as shear stress, drag, heat-transfer rate, and 
heat-transfer coefficient) associated with the free-convection 
flow. To this end, two parameters, a flow parameter and a 
heat-transfer parameter, are derived in appendixes C and D, 
respectively. 

The flow parameter 

r =^"(0) 

(4<7r x 3 )‘(*'„Mo/-X 2 ) 

is presented as a function of Prandtl number in figure 7. 
Thus, the various flow quantities for a given set of conditions 
can easily be computed by application of figure 7 . 

The local heat-transfer parameter 


x — J - t v-v 
(<?r x /4)‘ 

as determined here is given as a function of Prandtl number 
in figure 8. A calculation of the local Nusselt number from 
this equation for Pr =0.72 and 6r x ~ 10° yields a value of 
63.5, which indicates that large heat-transfer coefficients can 
also be obtained with free-convection flows. 



Fioube 6.— Comparison of largo plate experimental and theoretical velocity distributions 
for Prandtl number of 0.72. 


On the basis of a simplified theory (that is, by use of 
integrated momentum and energy equations and assumed 
velocity and temperature distributions), Eckert (see p. 162 
of ref. 7) obtained the approximate relation 

Nu _ 0.71 8(Pr)* 

(6-V x /4) i (0.952+Pr)* 
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Fioube 8.— Comparison of largo plato experimental and theoretical temperature distribu- 
tions for Prandtl number of 0.72. ^ 



Fioube 7.— Dimensionless flow parameter os function of Prandtl number. 
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Fioube 6.— Dim ensionl ess heat-transfer parameter as function of Prandtl number. 
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The curve representing this equation is also presented in 
figure 8, and it closely appro xim ates (to within about 10 
percent) the curve determined by the more exact considera- 
tions of this report over the entire Prandtl number range. A 
semiempixical equation as given in reference 11 relating the 
average (over the length A) Nusselt number to the Prandtl 
and Grashof numbers which has been used in the heat- 
transfer calculations up to the present is 


data from the oil experiments, though not so good, show 
reasonable agreement (maximum error of approximately 20 
percent) with the theoretical curve and good agreement, as 
is to be expected, with the semiempirical curve. The differ- 
ence between the theoretical values and the oil experiment 
results can possibly be due to the fact that the viscosity 
changes in oil are large even for small temperature differences 
or due to the end effects in the measurements. 


Nua,= 0.548 [(Pr)(?r]* 

The constant 0.548 pertains specifically to air; for oil it 
should be 0.555 (see ref. 12) and for mercury, approximately 
0.33 (see ref. 13). In order to obtain local values from the 
average ones given by the last equation, it is merely necessary 
to multiply the average values by 0.75. (The determination 
of this reduction factor of 0.75 is discussed in appendix D.) 
Thus in terms of the local quantities the semiempirical 
relation becomes 


or 


Nu= 0.411 [(Pr)6?r x ] * 


Nu 

«?r x /4) f 


0.581 (Pr)* 


The curve given by this equation is also presented in figure 
8 and the agreement with the theoretical curve is very good 
for Prandtl numbers near unity, not so good for large Prandtl 
numbers, and very poor for the small Prandtl numbers. Of 
course, changes of the constants in the semiempirical relation 
as previously discussed for the largo or small Prandtl number 
cases (oil and mercury, respectively, for example) would 
cause the semiempirical curve to approximate the theoretical 
curve more closely. The values of the heat-transfer param- 
eter obtained experimentally for mercury (Pr=0.03), air 
(Pr= 0.72), water (Pr= 7), and oil (Pr= 75.5, 115. 190, 224, 
275, 318, 368, and 442) are here reduced by the factor 0.75 
from the average values reported. The value for mercury 
is an average taken of four readings from a curve, since this 
experiment was the only one not reported in tabular form. 
From figure 8 it can be seen that all of the experimental 
values except those for the oil experiments are in very good 
agreement with the theoretically computed values. The 


CONCLUSIONS 

An analysis was made of the free-convection flow about a 
flat plate oriented in a direction parallel to that of the 
generating body force under the prime assumption that the 
relative temperature difference is small. It was found that 
the Grashof number was the principal factor determining 
the type of flow and that for large Grashof numbers the flow 
was of the boundary-layer type. The theoretical develop- 
ment was then continued to consider only the cases of largo 
Grashof number because these are of most importance in 
aeronautics. 

Velocity and temperature profiles for Prandtl numbers of 
0.01, 0.72, 0.733, 1, 2, 10, 100, and 1000 were computed on 
the basis of a constant body force and plate temperature 
and agreement with experiments where the fluid was air 
(Prandtl number of 0.72) was good. It was also demon- 
strated that velocities and Nusselt numbers of the ordor of 
magnitude of those obtained in forced-convection could be 
obtained in free-convection flows. 

A flow parameter and a heat-transfer parameter which are 
functions of the Prandtl number alone were derived. Calcu- 
lations of the important physical quantities such as shear 
stress, heat-transfer rate, and the like can be computed from 
these parameters. Values of the heat-transfer parameter 
obtained from an approximate theoretical development and 
from experiments compared with values computed from the 
present development showed good agreement over a wido 
range of Prandtl number (0.01 to 1000). It is shown that 
the commonly used semiempirical relation for the heat- 
transfer coefficient will yield good results only in restricted 
Prandtl number ranges. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, October 8, 1951 


APPENDIX A 


SYMBOLS 


The following notation is used in this report: 


A, in) ft An) 

C M , Di w 

F 

U 

y.v 


coefficients in numerical differentiation and 
integration formulas 
specific heat at constant pressure 
dimensionless velocity function 
components of body force per unit mass, 
i= 1, 2, 3 ■ 

negative of A-component of body Jorce per 
unit mass 


Or 

Gr x 

g 

H 

h 

K 


Grashof number, — “ — 5^ 

M=> 

Grashof number based on A 
gravitational force per unit mass (or accelera- 
tion due to gravity) 
dimensionless temperature function 
heat-transfer coefficient 
isothermal compressibility coefficient, 
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k 

l 

m, t 

N 

Nu 

Nu a , 

P 

Pr 

R 

r, s, t 

T 

U t 

u t 

u 

V 


Xt 
a 'i 
x 

Y 

V 


thermal-conductivity coefficient 
characteristic length 
arbitrary exponents 

a number, defined following equation (26) 
Nusselt number, hX/k 
average Nusselt number 
pressure 
Prandtl number 
gas constant 
arbitrary exponents 
absolute temperature 
velocity components, i= 1, 2, 3 
dimensionless velocity components, i= 1, 2, 3 
dimensionless velocity component in z-direc- 
tion 

dimensionless velocity component in y-direc- 
tion 

Cartesian coordinates, i= 1, 2, 3 
dimensionless Cartesian coordinates, i=l, 2, 3 
dimensionless Cartesian coordinate 
Cartesian coordinate 
dimensionless Cartesian coordinate 


P 

y 

A 

£ 

n 

e 

K 

M 

V 

p 

a 

T 

<P 

'P 

Subscripts: 


coefficient of volumetric expansion, 



ratio of specific heats 
Laplacian operator 

relative temperature difference, /3(T 0 — T m ) 

similarity variable 

dimensionless temperature function 

step size used in numerical calculations 

absolute viscosity 

kinematic viscosity 

density 

dimensionless pressure function 
shear stress 

dimensionless density function 
stream function 


i, j Cartesian tensor and summation subscripts 

8 denotes evaluation at static conditions (e=0) 

0 denotes evaluation at plate surface 

03 denotes evaluation at undisturbed conditions 

Subscript notation is used to denote partial differentiation. 
Superscripts: 

Primes denote ordinary differentiation. 

Bars (as cr or y) denote transformed dimensionless quantities. 


APPENDIX B 

NUMERICAL SOLUTION OF SIMPLIFIED BOUNDARY- VALUE PROBLEM 

By Lynn U. Aibbrs 


The method is presented herein by which solutions to the 
boundary-value problem 


F"'+3FF"—2F' 2 +H=0 

(Bl) 

H"+ZPrFH'=0 

<B2) 

F(0)=F'(0)=0 17(0) = 1 


^'(03) =#(<*> )=0 



wore obtained for the cases of Pr equal to 0.01, 0.72, 0.733, 
1, 2, 10, 100, and 1000. This discussion will enable the 
results to be clearly evaluated and will perhaps serve as a 
guide in the numerical solution of similar problems. 

Each of the cases of the problem has a solution for a 
particular set of values for F"(Q) and H'( 0), hereinafter 
called eigenvalues. The basic approach to the problem was 
to estimate the eigenvalues and to integrate out from zero, 
obtaining functions which satisfied equations (Bl) and (B2) 
at each step. The integration was continued until the func- 
tions F' and H behaved in a fashion inconsistent with the 
boundary values at infinity; for example, when they became 
negative or diverged to infinity. Improved estimates of the 
eigenvalues were then made on the basis of the results of 
preceding runs and the process was repeated successively 
until a solution was obtained. 


Modifications required to overcome specific obstacles will 
be discussed after sufficient details of the basic procedure 
have been given. Then an evaluation of the accuracy of the 
numerical results will be made. 

The integration process consists of two parts, a starting 
phase and an extension phase. The starting phase begins 
with an estimate of the eigenvalues F"{ 0) and H' (0) and a 
decision on the step size k to be used. It continues with an 
iterative process of alternately computing F"' and 27" at 
the first four points and integrating them by five-point 
formulas. This process and that in the extension phases 
are so arranged that the differential equations are satisfied 
at each integral multiple of the step size. 

The extension phase used preceding data to integrate step 
by step beyond the fourth point. Diagrams of both phases 
will be given after a few preliminary explanations. 

All integration formulas used are based on the same idea. 
If a function, for example, F'", is known at five points, there is 
a unique fourth-degree polynomial which agrees with it at 
these five points. Moreover, if the successive antideriva- 
tives (integrals) F", F', and F of F'" are known at one 
point, there are unique fifth-, sixth-, and seventh-degree 
polynomials which are successive, antiderivatives of this 
fourth-degree polynomial and which agree with F", F', and 
F, respectively, at the one point. It is then a simple algebra 
problem to deduce from the values of F"' at five points and 
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F, F' , and F" at a single point the values of any of these 
four polynomials at any point. These results will approxi- 
mate the functions F, F', F", and F'” to a degree dependent 
on step size, the relative positions of the points in question, 
and the magnitude of the fifth derivative of F"' in the 
neighborhood of these points. 

The preceding algebra problem can be presolved in all 
situations that arise in the starting and extension phases of 
the present problem and specific integration formulas may 
be deduced. These formulas are discussed in the next 
paragraph. 

Let F'" be denoted at five successive points by F 0 '", 
Fx'”, F a '", Fa'", and F"'. In the starting phase, these 
points are 0, k, 2k, 3k, and 4 k, and F 0 , F 0 ', and F 0 " are also 
known. Then the five sets of formulas required in the start- 
ing phase are . 



4 



A {J W F/" 


i=l, 2, 3, 4 (B3) 


H l '=H 0 '+ T ^±A t ^H/' 

]- o 

F i '=F 0 '+i K F 0 "+~ r pA {J <”F/" 
3 i =H 0 +iKE a '+ 1 ^T l A { ^H ) " 

J-Si jva 0 


i=l, 2, 3, 4 (B4) 
i= 1, 2, 3, 4 (B5) 


i= 1, 2, 3, 4 (B6) 


F i =F 0 +iKF 0 '+^~ F"- 




i=l, 2, 3, 4 (B7) 


where the superscripts on the A i] w and Z>< (B) refer to the 
order of integration. 

The constants A t] w and D ( w may be read from the fol- 
lowing tables: 

For A w (1) and D t ™: 


A«( 0 



■m 




■ 

D<V) 

■ 

m 

2 

3 

D 


1 

251 

m 

-264 

106 

EI 

720 

2 

29 

124 

24 

4 


90 

3 

27 

102 

72 

42 


80 

4 

14 

64 

24 

64 

14 

45 


For A,, 0 * and D t m : 


AfK 1 ) 


iX 






DiQO 

0 

1 

2 

3 

4 


1 

367 

540 

-282 

no 

-21 

1440 

2 

i 53 

144 

-30 

10 

El 

90 

3 

441 

1404 

162 

180 


480 

4 

56 

192 

48 

64 

1 

45 


For A i} G) and D t m : 




\ 1 






Dip) 

0 

1 

2 

3 

4 


1 

1017 

1070 

-618 

258 

-47 

10080 

2 

331 

664 

-240 

104 

-19 

630 

3 

1431 

3726 

-486 

460 

-81 

1120 

4 

744 

2176 

96 

384 

-40 

316 


It is now possible to diagram the steps of the starting phase 
of the integration. If each bar above a function denotes an 
improved estimate of it, and the first estimates of F\" , Fa" , 
Fa", and F"‘ are all equal to F 0 '", and similarly for the 
H", then the starting phase diagrams are 

(1) ( F 0 , Fo', F 0 ", Fa'", Fi", Fa'", F 3 ", F t '")->F U Fx', F x " 
(This diagram means that the values in parentheses are 
used ,with appropriate integration formulas from (B3) to 
(B7) to obtain F, F', and F” at i?=/c.) 

(2) (Ho, Ho', Ho", H", Ha", H", H<")^H U H t ' 

(3) (F x ; Ft', F x ", Hx, H'^Fx'", Hx" 

(The preceding diagram means that the values in parentheses 
are substituted in the differential equations (Bl) and (B2) 
to obtain F'" and H" at 17 —k.) 

(4) (F 0 , F 0 ', F 0 ", Fo'", lx'", Fa'", Fa'", F t "')-+Fa, Fa', F>" 

(5) (Ho, Ho', Ho", Hx", Ha", Ha", H")->H, Ha' 

(6) (Fa, Fa', Fa", H, H/)->Fa"' , Ha" 

(7) (Fo, Fo', Fo", Fo'", Fx'", la"', Fa'", F t '")->F 3 , F a ', F 3 " 

(8) (Ho, Ho', Ho", Hx", Ha", H", H")^H 3 , H 

(9) (Fa, Fa', Fa", H, H t ')-»Fa'", Ha" 

( 10 ) (Fo,Fo',Fo",Fo''',lx"',!t'",l,"',Fr)->F i ,F i ',F i " 

(11) (Ho, Ho', Ho", Hx", Ha", Ha", H' 

(12) (Ft, Ft', Ft", Ht, Ft')->1 "', Ht" 

It may be noted here that all four values of F'" and H" 
have been improved, and further improvement will require 
iteration of steps 1 to 12. The start of the second iteration 
is diagrammed as follows: 

(13) (Fo,Fo',Fo'',Fo''',!x'",!a'",l 3 '",Ht'")->lx,!i',li" 

(14) (H, H', H", Hx", Ha" JIz" ,Jdt")~>Hx, Hx' 

(15) (Fx, Fx', Fx", Hx, Hx')^Fx'", Hx" 

(16) (Fo,Fo',Fo",Fo"',Fx'",Fa'",Fa'",lt'")^Fa,lV,la" 


Successive sets of 12 steps are performed until the values of 
Fx'” and H " no longer change. 

On the IBM Card-Programmed Electronic Calculator, a 
deck of punched cards 2 inches thick sufficed to perform steps 
1 to 12. Three runs of this starter deck at 3 minutes per 
run accomplished complete convergence in most cases. At 
the end of the starting process there have been computed and 
stored Ft, Ft, F" , Ht, and S/, and final estimates of Fx"', 
Fa'", Fa'", Ft'", Hx", Ha", H", and H". 





























FREE-CONVECTION PLOW AND HEAT TRANSFER ABOUT A FLAT PLATE 75 


The extension phase has now been reached. It used a 
different set of integration formulas based on the same gen- 
eral ideas as equations (B3) to (B7). If Fo" , Ft" , Fo" , 
Ft'", and Ft" now designate F'" at any five successive 
points, and the subscript 5 denotes the next point, 

Ft"=Ft"+-^p o B^F/" 

(B8) 

Ht'=Ht'+-^p o B/»H ] " 

<B9) 

Ft'=Ft'+KFt"+-^ g B^F/" 

(B10) 

Ht=Ht+ kH 7 g B,® H" 

(Bll) 

F s =Ft+KFt'+£ Ft"+-~±B^F'" 

(B12) 


where the B/ n) and <7 tn) are given in the following table: 


B/W 

o> 


0 

1 

2 

3 

4 

1 

251 

—1274 

2616 

-2774 

1901 

720 

2 

136 

-692 

1449 

-1596 

1427 

1440 

3 

<10 

-2116 

4476 

-50S4 

6674 

20160 


The extension phase may then be diagrammed simply as 
follows: 

(1) (F<, F/, Ft", Fo"', Ft'", F*"', F/", P/")->P c , F t ', F 6 " 

(2) (Bt, Ht', Ho", Ht", Ho", Ho", fl/')-*#., HJ 

(3) (F b , FJ, Fo", Hs)->F b '", Ht" 

The values of the functions at the next point are computed 
in similar manner, where the latest sets of five values of 
F'" and H" are used. This process advances step by step 
toward infinity. 

The extended deck of punched cards was about 3 inches 
thick and took a little over 3 minutes per run. For Pr=0.72, 
a step size of 0.1 was used, the starting phase took 10 minutes, 
and the extension phase, about 30 minutes. When it is 
realized that about 11,000 operations were performed in the 
40 minutes per run, it may be seen that solution of the 
present problem would have been prohibitively difficult on 
desk-type calculators. Simplifications in method would have 
sacrificed accuracy or required smaller step size. 

In two-point boundary-value problems where one point is 
infinity, some problems of judgment are involved as to where 
infinity is, and as to when a satisfactory appro xima tion to a 
solution has been obtained. IiL most cases this question was 
settled for the present problem by c alling a ran satisfactory 
when it fell between two runs for which F' and H did not 
differ at important points in the fourth decimal place, and 


for which F' and H flattened out at zero, correct to four 
decimal places. 

Certain difficulties were met in the attempt to use the 
basic procedure previously discussed. These necessitated 
certain modifications. 

For Pr= 2, 10, 100, and 1000, H would settle down to 
zero at an early stage; but while F' was still coming down, 
H" would begin to oscillate and these oscillations increased 
and fed back into all other functions. This trouble was 
avoided by the following modifications: It is a consequence 

of equation (B2) that 

H'( V )=H'( 0) exp (-3Prj^F(t)dl) 

=H'(i) — k) exp 3PrJ* F(t)dtJ (B13) 

The extension phase was modified to require the additional 
integration formulas 

£/ S AiFi (B14) 

Ht=Ht+^^AiH t ' (B16) 

where Ai=— 19, A 2 =106, yl 3 =— 264, AU=646, and ^=261. 

These formulas were used along with equations (B8) to 
(B10) according to the following diagram: 

( 1 ) (F h Ft', Ft", Fo'", Ft'", F 2 '", Fo"', FJ")-+F S , F s ', Ft" 

(2) (Fi, F 3 , F z , Ft, Ft, Ht')->H' by means of (B14) and 
(B13) 

(3) (H/, H/, Ho', Ht', Ht', Ht)->H 5 by means of (B15) 

( 4 ) (Ft, Ft', Ft", Ht)->Ft'" 

The value F 0 '" is discarded and F"' at the last five points 
is used to repeat the whole process again and again ad 
infinitum. As long as F stays positive, H' is guaranteed to 
approach zero and H will flatten out to some value and not 
oscillate. 

For Pr=0.01, 0.72, 0.733, $nd 1, the F"' began to oscillate 
at an advanced point and these oscillations grew and fed into 
the other functions. For all cases but Pr=0.01, the oscilla- 
tions appeared very late, near the end of the run, and a 
suitable halving of step size when oscillation was detected 
in the fourth differences of F'" was sufficient to avoid the 
difficulty. But for the 0.01 case, oscillations of F'" appeared 
early in the run, namely, soon after the peak in F. These 
oscillations were found to be step-size connected, so that 
reduction of the step to 0.02 avoided them. Even then 
oscillations in F r " would begin to appear every 25 steps or 
so, and these were smoothed out regularly by repeated runs 
of a deck similar to the starter deck. Each run under these 
conditions took about 16 hours, making this the most diffi- 
cult case to solve. 
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APPENDIX C 

DERIVATION OF FLOW PARAMETER 


By definition the shear stress is given by 


t=m 0 



(Cl) 


Substitution of this expression into equation (Cl) yields the 
flow parameter 


=^"( 0 ) 

(±Gr x *)\ v ^IX*) 


To express (bUfbY) 0 in terms of the known function F{n), 
use can he made of equations (60) and (62). Then 


g-cWji 


Note that from the general derivation, the flow parameter 
contains the viscosity evaluated at two different points. 
Becall, however, that the analysis has shown that to a first 
approximation the variation of viscosity with temperature 
can be neglected. Thus the viscosity can be taken as 
constant in the entire flow field. 


APPENDIX D 

DERIVATION OF HEAT-TRANSFER PARAMETER 


The local Nusselt number is defined as 



-X 

(To-TJ 



(Dl) 


To express (bT/bY) 0 in terms of the known function Hirf), 
use is made of equations (61) and (62). Thus 


bT (T 0 -T J 
bY~ X 



The heat-transfer parameter as given by equation (D2) 
is, as was previously stated, a local parameter. It is often 
desired to compute the average (over the length X) value 
of this parameter. To this end, the Nusselt number (as 
given in equation (Dl)) must he defined in terms of an 
average heat-transfer coefficient and the quantity thus ob- 
tained must then be integrated over the length X and 
divided by X. This procedure yields the result 


Substitution of this expression into equation (Dl) yields 
the heat-transfer parameter 


Nu 

(Gr x / 4) f 


S'i 0 ) 


(D2) 


(Nu)„ 

It is from this last equation that the 0.76 reduction factor 
previously discussed was obtained. 
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TABLE I.— FUNCTIONS F AND H AND DERIVATIVES FOR VARIOUS PRANDTL NUMBERS 


(a) Prandtl number, 0.01 


V 

F 

F 

F " 

H 

H ' 

V 

F 

F * 

F " 

H 

H' 

0 
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0.0000 

0.9862 

1.0000 
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.9919 

—.0812 
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.9838 

—.0812 

4.20 

L 9480 

.4012 

—.0612 

.6712 

—.0721 

.3 

.0399 

.2516 

.6942 
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2.1780 

.3660 

—.0564 

.6287 

—.0695 

1.2 

.4480 

.5631 

.0839 

.9028 
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1.1703 
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.7227 

—.0749 

11.80 

a 7086 
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(b) Prandtl number, 0.T2 (c) Prandtl number, 0.733 
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TABLE I.— FUNCTIONS F AND H AND DERIVATIVES FOR VARIOUS PRANDTL NUMBERS— Continued 


(d) Prandtl number, 1 
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-.6371 

.7 

.1001 

.2406 

.1008 

.0116 

—.5221 

.8 

.1308 

.2481 

.0516 

.5602 

-.5038 

.9 

.1558 

.2511 

.0093 

.6109 

-.4828 

1.0 

.3809 

.2502 

-.0263 

.4638 

-.4589 

LI 

.2058 

.2401 

-.0557 

.4192 

-.4330 

L 2 

.2300 

.2393 

-.0793 

.3772 

-.4050 

1.3 

.2635 

.2301 

-.0975 

.3381 

-.3772 

L 4 

.2701 

.2199 

-. 1110 

.3018 

—.3484 

1.5 

.2976 

.2083 

-.1203 

.2884 

-.3197 

1.6 

.3177 . 

.1900 

-.1280 

.2379 

-.2915 

L 7 

.3367 

.1832 

-.1287 

.2101 

-.2842 

L 8 

.3643 

.1708 

-.1288 

.1860 

— . 23 S 2 

1.9 

.3707 

.1575 

-.1288 

.1624 

-.2136 

2.0 

.3859 

.1450 

-.1233 

.1422 

-.1907 

2.1 

.3997 

.1239 

-.1185 

.1242 

-.1695 

2.2 

.4125 

.1213 

-. 1127 

.1082 

—.1601 

2.3 

.4240 

.1104 

-.1064 

.0941 

— .1324 

2.4 

.4340 

.1001 

-.0997 

.0817 

-.1104 

2.5 

.4441 

.0901 

-.0928 

.0708 

-. 1020 

2.6 

.4527 

.0815 

-.0859 

.0013 

-.0892 

2.7 

.4604 

.0733 

-.0791 

.0529 

—.0777 

2.8 

.4673 

.0057 

-.0725 

.0457 

-.0676 

2.9 

.4736 

. 05 S 8 

-.0062 

.0392 

-.0537 

3.0 

.4791 

.0524 

-.0602 

.0339 

-.0509 

3.1 

.4841 

.0107 

-.0546 

- .0291 

- 0441 

3.2 

.4885 

.0116 

-.6493 

.0250 

-.0381 

3.3 

.4924 

.0368 

-.0444 

.0215 

-.0329 

3.4 

.4959 

.0328 

-.0399 

.0185 

-.0283 

3.6 

.5010 

.0254 

-.0321 

.0130 

—.0210 

3.8 

.5061 

.0197 

-.0255 

.0009 

-.0155 

4.0 

.5090 

.0151 

-.0202 

.0072 

-.0115 

4.2 

.5122 

.0110 

-.0158 

.0053 

-.0084 

4.4 

.6143 

.0087 

-.0124 

. .0038 

-.0002 

4.6 

.5158 

.0060 

-.0090 

.0027 

-.0046 

4.8 

.5169 

.0019 

-.0075 

.0020 

-.0034 

5.0 

.5177 

.0035 

-.0057 

.0014 

-.0024 

5.2 

.5183 

.0025 

-.0014 

.0010 

-.0018 

6.5 

.5189 

.0014 

-.0029 

.0006 

-.0011 

6.0 

■ 51 M 

.0001 

-.0014 

.0002 

-.0005 

6.25 

.6194 

.0000 

-.0010 

.0000 

-.0004 


(1) Prandtl number, 10 


(e) Prandtl number, 2 


V 

F 

F 

F" 

H 


0 

0.0000 

0.0000 

a 5713 

1.0000 

-a 7165 

.1 

.0027 

.0523 

.4749 

.9284 


.2 

.0101 

.0952 

.3861 

.8569 


.3 

.0215 

.1297 

.3049 

.7858 


.4 

.0358 

-1665 

.2318 

.7157 

-.0949 

.5 

.0525 

.1763 

.1666 

.6470 

-.6763 

.6 

.0709 

.1901 

.1096 

.5805 

-.0523 

.7 

.0961 

.1965 

.0602 

.6168 

-.6215 

.8 

.1104 

.2024 

.0185 

.4564 

sag 

.9 

.1307 

.2024 

-.0101 

.3999 

W* 

L 0 

>.1608 

.1994 

-.0440 

.3476 

’■SI A i l»« 

LI 

.1705 

.1938 

-•0659 

.2999 


1.2 

.1896 

.1861 

—.0824 

.2568 

.WrrM 

L 3 

.2077 

.1775 

-.0942 

.2183 

-.3619 

L 4 

.2250 

.1677 

-.1020 

.1844 

-.3178 

L 5 

.2412 

.1672 

-.1063 

.1547 

-.2763 

L 6 

.2561 

.1405 

-.1078 

.1290 


L 7 

.2706 

.1357 

-.1071 

.1069 


1.8 

.2836 

.1252 

-.1040 

.0882 

-.1721 | 

L 9 

.2956 

.1149 

-.1008 

.0724 

-.1440 

2.0 

.3066 

.1050 

-.0960 

.0592 

Ei r 

2.1 

.3166 

.0957 

-.0906 

.0482 

l| J ll 

2.2 

.3257 

.0869 

-.0849 

.0391 

i :W; 

2.3 

.3340 

.0787 

-.0789 

.0316 

IB® i l ;Vf 

2.4 

.3416 

.0711 

-.0729 

.0254 

iMJy 

2.5 

.3483 

.0611 

-.0071 

.0204 

|M/I 

2.6 

.3543 

.0677 

-.0614 

.0104 


2.7 

.3598 

.0518 

-.0560 

.0131 

iW« 

2.8 

.3647 

.0465 

-.0609 

.0105 


2.9 

.3691 

.0417 

-.0461 

.0083 

-.0100 

3.0 

.3731 

.0373 

-.0416 

.0066 

-.0152 

3.1 

.3760 

.0333 

-.0375 

.0053 

-.0121 

3.2 

.3798 

.0298 

-.0338 

.0042 

-.0097 

3.3 

.3826 

.0266 

-.0303 

.0033 

-.0077 

3.4 

.3851 

.0237 

-.0272 

.0020 

-.0061 

3.0 

.3893 

.0188 

-.0218 

.0017 

-.0038 

3.8 

.3927 

.0149 

-.0174 

.0010 

-.0024 

4.0 

.3953 

.0018 

-.0138 

.0007 

-.0016 

4.2 

.3974 

.0094 

-.0109 

.0004 

-.0009 

4.4 

.3991 

.0074 

-.0086 

.0003 

y&ra I’n'i 

4.6 

.4004 

.0059 

-.0068 

.0002 


4.8 

.4016 

.0047 

-.0054 

.0001 

iVi/ 

5.0 

.4023 

.0037 

-.0042 

.0001 

H m y 

5.4 

.4035 

.0024 

-.0020 

.oooi 

ijM 

5.8 

.4043 

.0015 

-.0016 

.0000 

^B t**u 

6.4 

.4049 

.0008 

-.0008 

.0000 

jilt- iffl 

7.0 

.4053 

.0005 

-.0001 

.0000 

f^R l,l*H 

8.0 

.4056 

.0002 

-.0001 

.0000 

iVi’i 

9.0 

.4058 

.0001 

.0000 

.0000 

l^R iVn 

10.0 

.4050 

.0001 

.0000 

.0000 

ijil 

11.0 

.4059 

..0001 

.0000 

.0000 

.0000 


(g) . Prandtl number, 100 


V 

F 

F 

F ‘ 

H 

H ' 


V 

F 

F 

F ’ 

H 

H ’ 

0 

0.0000 

a oooo 

0.4192 

L 0000 

— L 1694 


0.000 

0.0000 

0.0000 

0.2517 

L 0000 

- 2. 191 

.1 

.0019 

.0371 

.3251 

.8832 

— L 1671 


.025 

.0001 

.0060 

.2274 

.9452 

-2. 101 

.2 

.0071 

.0654 

.2428 

.7670 

- L 1521 


.050 

. .0003 

.0114 

.2044 

.8905 

- 2. 188 

.3 

.0147 

.0801 

.1723 

.6534 

- L 1155 


.075 

.0006 

.0162 

.1828 

.8369 

- 2.180 

.4 

.0241 

.1003 

.1134 

.5448 

- L 0526 


.100 

.0011 

.0205 

.1620 

.7815 

- 2. 106 

.5 

.0340 

.1091 

.0655 

.4437 

-.9040 


• .125 

.0017 

.0244 

.1438 

.7270 

- 2.144 

.0 

.0458 

.1137 

.0279 

.3527 

—.8545 


.160 

.0023 

.0277 

.1283 

.6744 

— Z 113 

.7 

.0573 

.1150 - 

-.0011 

.2733 

-.7322 


.175 

.0031 

.0307 

.1101 

.6221 

- 2.071 

.8 

. 06 S 7 

.1137 

-.0221 

.2064 

-.0061 


.200 

.0039 

.0332 

.0952 

.5709 

- Z 018 

.9 

.0800 

.1107 

-.0307 

.1519 

-.4849 


.225 

.0047 

.0355 

.0816 

.5213 

— L 954 

LO 

.0908 

.1066 

-.0462 

.1090 

-.3753 


.250 

.0050 

.0373 

.0692 

.4733 

- L 8 S 0 

LI 

.1012 

.1016 

-.0518 

.0763 

-.2813 


.275 

.0006 

.0389 

.0580 

.4273 

— L 780 

L 2 

.1111 

.0963 

-.0545 

. 0522 ’ 

-.2045 


.300 

.0070 

.0402 

.0479 

.3836 

- 1.704 

L 3 

.1205 

.0908 

-.0552 

.0349 

-. 1445 


.325 

.0080 

.0413 

.0389 

.3422 

- 1.604 

L 4 

.1293 

.0853 

— .0544 

.0228 

-.0993 


.350 

.0096 

.0422 

.0309 

.3034 

- L 498 

L 6 

.1370 

.0799 

-.0527 

.0146 

-.0665 


.375 

.0107 

.0429 

.0238 

.2874 

- L 38 S 

L 6 

.1453 

.0748 

-.0505 

.0092 

-.0435 


.'400 

.0118 

.0434 

.0176 

.2341 

- L 270 

L 7 

.1525 

.0699 

-.0480 

.0056 

-.0278 


.425 

.0129 

.0438 

.0123 

.2030 

— L 163 

L 8 

.1593 

.0652 

—.0453 

.0034 

-.0174 


.450 

.0140 

.0440 

.0078 

.1769 

- L 052 

L 9 

.1656 

.0608 

-.0427 

.0020 

-.0107 


.475 

.0151 

.0442 

.0030 

.1610 

-.9130 

2.0 

.1714 

.0567 

-.0400 

.0012 

-.0065 


.500 

.0162 

.0442 

.0002 

.1287 

-.8393 

2.1 

.1769 

.0538 

-.0375 

.0007 

-.0038 


.525 

.0173 

.0442 

-.0026 

.1089 

-.7404 

2.2 

.1820 

.0491 

-.0351 

.0004 

-.0022 


.560 

.0184 

.0441 

-.0060 

.0910 

-.6178 

2.3 

.1808 

. 045 S 

-.0328 

.0002 

-.0013 


.575 

.0195 

.0439 

-.0070 

.0706 

-. 6021 

2.4 

.1912 

.0426 

-.0306 

.0001 

-.0007 


.600 

.0200 

.0437 

-.0087 

.0634 

-.4837 

2.5 

.1953 

.0396 

-.0286 

.0001 

-.0004 


.65 

.0227 

.0432 

-.0111 

.0427 

-.3490 

2.6 

.1991 

.0369 

-.0287 

.0000 

-.0002 


.70 

.0249 

.0426 

-.0128 

.0280 

3440 

2.7 

.2027 

.0343 

-.0249 

.0000 ■ 

-.0001 


. 76 . 

.0270 

.0420 

-.0135 

.0178 

-. 1657 

2.8 

.2060 

.0319 

-.0232 

.0000 

-.0001 


.80 

.0291 

.0413 

-.0140 

.0111 


2.9 

.2090 

.0297 

-.0216 

.0000 

.0000 


.85 

.0311 

.0406 

-.0142 

.0067 


3.0 

.2119 

.0270 

-.0201 

.0000 

.0000 


.90 

.0331 

.0399 

-.0142 

.0039 


3.1 

.2140 

.0256 

-.0187 

.0000 

.0000 


.96 

.0351 

.0392 

-.0141 

.0022 


3.2 

.2170 

.0238 

-.0174 

.0000 

.0000 


LOO 

.0371 

.0385 

-.0140 

.0012 

-. 0149 

3.3 

.2193 

.0221 

-.0162 

.(fooo 

.0000 


L 10 

.0408 

.0371 

-.0130 

.0004 


3.4 

.2215 

.0206 

-.0151 

.0000 

.0000 


L 20 

.0445 

.0357 

-. 0132 

.0001 

[JIB i*i) ^ 

3.6 

.2253 

.0178 

0131 

.0000 

.0000 


LSO 

.0480 

.0344 

-.0128 

.0000 

l§B ijjV 

3.8 

.2286 

.0153 

-.0113 

.0000 

.0000 


L 40 

.0614 

.0332 

-.0124 

.0000 

iVn 

4.0 

.2314 

.0132 

-.0098 

.0000 

.0000 


LSO 

.0540 

.0320 

-.one 

.0000 


4.2 

.2339 

.0114 

-.0084 

.0000 

.0000 


L 00 

.0578 

.0308 

-.0116 

.0000 

||B M . l,l 

4.4 

.2360 

.0098 

-.0073 

.0000 

.0000 


L 70 

.0608 

.0297 

-.0112 

.0000 

f|B iVl l 

4.6 

.2379 

.0085 

-.0063 

.0000 

.0000 


LSO 

.0637 

.0286 

-.0108 

.0000 

f|B 

4.8 

.2394 

.0073 

-.0054 

.0000 

.0000 


LOO 

.0065 

.0275 

-.0104 

.0000 

i|B i|iVi 

5.0 

.2408 

.0063 

-.0047 

.0000 

.0000 


2.0 

.0692 

.0285 

-.0101 

.0000 

gB iViti 

5.4 

.2430 

.0047 

-.0035 

.0000 

.0000 


21 

.0718 

.0255 

-.0097 

.0000 

gB nVi 

5.8 

.2446 

.0035 

-.0028 

.0000 

.0000 


22 

.0743 

.0245 

—.0094 

.0000 

f|B iVfl 

8.2 

.2458 

.0026 

-.0019 

.0000 

.0000 


23 

.0767 

.0230 

-.0091 

.0000 

|JB nVi 

7.0 

.2474 

.0014 

-.0011 

.0000 

.0000 


24 

.0790 

.0227 

-.0088 

.0000 

IB nVi 

8.0 

.3484 

.0007 

-.0005 

.0000 

.0000 


20 

.0834 

.0210 

-.0082 

.0000 

||B iVili 

9.0 

.3489 

.0003 

-.0002 

.0000 

.0000 


28 

.0874 

.0105 

-.0070 

.0000 

M rriti 

10.0 

.3491 

.0002 

-.0001 

.0000 

.0000 


3.0 

.0912 

.0180 

-.0071 

.0000 

.0000 
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TABLE I.— FUNCTIONS F AND H AND DERIVATIVES FOR VARIOUS PRANDTL NUMBERS— Concluded 

(g) Prandtl number, 100— Concluded (h) Prandtl number, 1000 


V 

F 

F* 

F" 

H 

R' 


V 

F 

F> 

F" 

H 

H 7 

3.2 

0.0917 

a 0166 

- 0.0066 

0.0000 

a oooo 


a 

a oooo 

a oooo 

a 1450 

1.0000 

- 3.966 

3,4 

.0979 

.0151 

-.0001 

.0000 

.0000 


.025 

.0000 

.0033 

.1212 

.9009 

- 3.962 

3.6 

.1008 

.0142 

-.0056 

.0000 

.0000 


.050 

.0002 

.0061 

IIK Liu! 

.8021 

- 3.933 

3.8 

.1035 

.0131 

-.0052 

.0000 

.0000 


.075 

.0003 

.0083 

iJB Si 

.7010 

- 3. 801 

4.0 

.1001 

.0121 

-. 0019 

.0000 

.0000 


.100 

.0006 

.0102 


.6090 

- 3.731 

4.4 

.1105 

.0103 

-. 0042 

.0000 

.0000 


.125 

.0008 

.0116 

hb iVi*. 

.5186 

- 3.638 

4.8 

.1143 

.0038 

-.0036 

.0000 

.0000 


.160 

.0012 

.0127 

|||B i?:f 

.4332 

- 3.283 

6.2 

.1176 

.0074 

-.0031 

.0000 

.0000 


.175 

.0015 

.0136 

i||B iffii 

.3549 

- 2.975 

6.6 

.1203 

.0063 

-.0026 

.0000 

.0000 


.200 

.0018 

.0142 

ilfB iVTiV 

.2847 

- 2.628 

6.0 

.1226 

.0053 

-.0022 

.0000 

.0000 


.225 

.0022 

.0140 

||B iTrfi 

.2236 

- 2.201 

6.6 

.1264 

.0011 

-.0018 

.0000 

.0000 


.250 

.0020 

.0149 

HB 

.1717 

- 1.893 

7.2 

.1270 

.0032 

-.0014 

.0000 

.0000 


.275 

.0029 

.0151 

||B mV' 

.1283 

- 1.641 

ao 

.1297 

.0022 

-.0010 

.0000 

.0000 


.300 

.0033 

.0152 

IB I ’ KV . 1 

.0944 

- 1.220 

9.0 

.1316 

.0014 

-.0007 

.0000 

.0000 


.325 

.0037 

.0153 

'JIB i*n ^ 

.0675 

-.9331 

10,0 

.1326 

■ OOOS 

-.0005 

.0000 

.0000 


.350 

.0041 

.0153 

IB uY; 

.0471 

-.7012 

11.0 

.1332 

.0001 

-.0003 

.0000 

.0000 


.375 

.0045 

.0162 

llfi riTP 

.0321 

-.5003 

12.0 

.1335 

.0002 

-.0002 

.0000 

.0000 


.400 

. 004 S 

.0152 

gsfi ffTT 


-.3696 

13.0 

.1336 

.0000 

-.0001 

.0000 

.0000 


.425 

.0052 

.0162 

3SB iTx 


-.2407 








.460 

.0066 

.0151 

hb TOi 


-.1645 








.476 

.0060 

.0150 

|e rj 2 


-.1066 








.500 

.0063 

.0150 

fIB ffi ?: 


-.0872 








.525 

.0067 

.0149 



-.0413 








.560 

.0071 

.0148 

HjE rff ? 


-.0245 








.575 

.0076 

. .0147 



-.0142 








.600 

.0078 

.0147 

IB 22 


-.0080 








.625 

.0082 

.0146 

||E y 2 


-.0044 








.800 

.0107 

.0141 

jpg J 2 


. .0000 








1.000 

.0136 

.0136 

gg mV / 


.0000 








1.40 

.0187 

.0125 

|sB $2 


.0000 








1.80 

.0235 

. .0116 

mS u/V 


.0000 




* 




2.20 

.0279 

.0108 

-.0022 

yHvmjV 

.0000 








2.60 

.0320 

.0098 

—.0020 


.0000 








3.0 

.0368 

.0090 

-.0019 


.0000 








3.6 

.0109 

.0080 

-. 0017 


.0000 








4.2 

.0464 

.0070 

-.0015 


.0000 








6.0 

.0505 

.0060 

-. 0012 


.0000 








5.8 

.0619 

.0050 

-.0011 


.0000 








7.0 

.0603 

.0039 

-.0008 


.0000 








8.0 

.0638 

.0032 

-.0007 


.0000 








iao 

.0691 

.0022 

-.0001 


.0000 








12.0 

.0727 

.0015 

-.0003 


.0000 








14.0 

.0752 

.0011 

-.0002 


.0000 








iao 

.0771 

.0008 

-.0001 

.0000 

.0000 








18.0 

.0788 

.0007 

.0000 

.0000 

.0000 








2a o 

.0798 

.0008 

.0000 

.0000 

.0000 








22.0 

.0809 

.0005 

.0000 

.0000 

.0000 


• 





23.6 

.0816 

.0005 

.0000 

.0000 

.0000 


o 





